SYSTEMS AND METHODS FOR DETERMINING SPECTRA USING DYNAMIC 
LEAST SQUARES ALGORITHMS WITH MEASUREMENTS FROM LED 
COLOR SENSOR 
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BACKGROUND OF THE INVENTION 

1. Field of Invention 

[0002] This invention relates to determining spectra based on non-spectral 

inputs. 

2. Description of Related Art 

[0003] Automatic on-line color calibration systems can be much more 
effective with an on-line color measurement system where a spectrophotometer may 
be mounted in the paper path of the moving copy sheets in the printer, preferably in 
the output path after fusing or drying, without having to otherwise modify the printer, 
or interfere with or interrupt normal printing, or the movement of the printed sheets in 
said paper path, and yet provide accurate color measurements of test color patches 
printed on the moving sheets as they pass the spectrophotometer. That enables a 
complete closed loop color control of a printer. 

[0004] A typical spectrophotometer gives color information in terms of 
measured reflectances or transmittances of light, at the different wavelengths of light, 
from the test surface. This spectrophotometer desirably provides distinct electric 
signals corresponding to the different levels of reflected light received from the 
respective different illumination wavelength ranges or channels. 
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[0005] Known devices capable of providing distinct electric signals 
corresponding to the different levels of reflected light received from the respective 
different illumination wavelength ranges or channels include a grating-based 
spectrophotometer made by Ocean Optics Inc., LED based sensors marketed by 
5 "ColorSawy" or Accuracy Microsensor; and other spectrophotometers by Gretag 
MacBeth (Viptronic), ExColor, and X-Rite (DTP41). However, those devices are 
believed to have significant cost, measurement time, target displacement errors, and/or 
other difficulties, for use in real-time printer on-line measurements. 

[0006] As used herein, unless otherwise specifically indicated, the term 

1 0 "spectrophotometer" may encompass a spectrophotometer, colorimeter, and 

densitometer, as broadly defined herein. The definition or use of such above terms 
may vary or differ among various scientists and engineers. However, the following is 
an attempt to provide some simplified clarifications relating and distinguishing the 
respective terms "spectrophotometer," "colorimeter," and "densitometer," as they may 

1 5 be used in the specific context of specification examples of providing components for 
an on-line color printer color correction system, but not necessarily as claim 
limitations. 

[0007] A typical "spectrophotometer" measures the reflectance of an 
illuminated object of interest over many light wavelengths. Typical prior 

20 spectrophotometers in this context use 1 6 or 32 channels measuring from 380 nm to 
730 nm or so, to cover the humanly visible color spectra or wavelength range. A 
typical spectrophotometer gives color information in terms of measured reflectances 
or transmittances of light, at the different wavelengths of light, from the test surface. 
(This is to measure more closely to what the human eye would see as a combined 

25 image of a broad white light spectra image reflectance, but the spectrophotometer 
desirably provides distinct electrical signals corresponding to the different levels of 
reflected light from the respective different illumination wavelength ranges or 
channels.) 

[0008] A "colorimeter" normally has three illumination channels, red, green 
30 and blue. That is, generally, a "colorimeter" provides its three (red, green and blue or 
"RGB") values as read by a light sensor or detector receiving reflected light from a 
color test surface sequentially illuminated with red, green and blue illuminators, such 
as three different color LEDs or one white light lamp with three different color filters. 
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It may thus be considered different from, or a limited special case of, a 
"spectrophotometer," in that it provides output color information in the trichromatic 
quantity known as RGB. 

[0009] Trichromatic quantities may be used for representing color in three 
5 coordinate space through some type of transformation. Other RGB conversions to 
"device independent color space" (i.e., RGB converted to conventional L*a*b*) 
typically use a color conversion transformation equation or a "lookup table" system in 
a known manner. 

[0010] A "densitometer" typically has only a single channel, and simply 

1 0 measures the amplitude of light reflectivity from the test surface, such as a developed 
toner test patch on a photoreceptor, at a selected angle over a range of wavelengths, 
which may be wide or narrow. A single illumination source, such as an IR LED, a 
visible LED, or an incandescent lamp, may be used. The output of the densitometer 
detector is programmed to give the optical density of the sample. A densitometer of 

1 5 this type is basically "color blind." For example, a cyan test patch and magenta test 
patch could have the same optical densities as seen by the densitometer, but, of 
course, exhibit different colors. 

SUMMARY OF THE INVENTION 
[0011] A multiple LED reflectance spectrophotometer, as in the examples of 

20 the embodiments herein, may be considered to belong to a special class of 

spectrophotometers which normally illuminate the target with narrow band or 
monochromatic light. Others, with wide band illumination sources, can be flashed 
Xenon lamp spectrophotometers, or incandescent lamp spectrophotometers. A 
spectrophotometer is normally programmed to give more detailed reflectance values 

25 by using more than 3 channel measurements (for example, 10 or more channel 
measurements), with conversion algorithms. That is in contrast to normal three 
channel colorimeters, which cannot give accurate, human eye related, reflectance 
spectra measurements, because they have insufficient measurements for that (only 3 
measurements). 

30 [0012] It is desirable for a printer color control system to dynamically 

measure the color of test patches on the printed output media "on line", that is, while 
the media is still in the sheet transport or paper path of a print engine, for real-time 
and fully automatic printer color correction applications. 



4 

[0013] For a low cost implementation of the color sensor, a multiple 
illuminant device is used as the illumination source, and has, for example, 8, 10, 12 
or 1 6 LEDs. Each LED is selected to have a narrow band response curve in the 
spectral space. Therefore, for example, ten LEDs would correspond to ten 
5 measurements in the reflectance curve. The LEDs, or other multiple illuminant based 
color sensor equivalent, e.g., lasers, are switched on one at a time as, for example, the 
measured media is passed through a transport of a printer. The reflected light is then 
detected by a photodetector and the corresponding voltage integrated and normalized 
with a white tile. 

1 0 [0014] To obtain a smooth curve similar to that of a Gretag 

spectrophotometer, linear or cubic spline algorithms could be used, which blindly 
interpolates the data points without the knowledge of the color space. Unfortunately, 
due to lack of measurements at wavelengths below 430nm and above 660nm (due to 
lack of LEDs at these wavelengths), extrapolation with 10 measurements can lead to 

1 5 errors. 

[0015] The systems and methods of this invention use the integrated sensor 
measurements to determine a fully populated reflectance spectra with reflectance 
values at specific wavelengths, even though some of the light sources may not 
produce spectral content at the distant ends of the visible spectrum. By using a 

20 reconstruction algorithm, based on the spectral characteristics of the illumination 

source and the color sensing system, the integrated multiple illuminant measurements 
from a non-fully illuminant populated color sensor are converted into a fully 
populated spectral curve. 

[0016] Algorithms according to this invention utilize a reference database that 

25 contains training samples that indicate reflectance spectra and their corresponding 
LED sensor output. A dynamic, Least Squares-based (DLS) spectral reconstruction 
algorithm is used to reconstruct spectra. The algorithm is "dynamic" because it gives 
greater importance to the data from the training samples in the neighborhood of the 
color sample under measurement. This is done using linear operators. 

30 [0017] These and other objects, advantages and salient features of the 

invention are described in or apparent from the following description of exemplary 
embodiments. 



BRIEF DESCRIPTION OF THE DRAWINGS 
[0018] Exemplary embodiments of the invention will be described with 
reference to the drawings, wherein like numerals represent like parts, and wherein: 

Fig. 1 is a functional block diagram illustrating an exemplary embodiment of a 
5 coloring system according to the invention; 

Fig. 2 is a flowchart illustrating a first exemplary method of determining 
spectra according to this invention; 

Fig. 3 is a flowchart illustrating a second exemplary method of determining 
spectra according to this invention; and 
10 Fig. 4 is a functional block diagram illustrating an exemplary embodiment of a 

color detection system according to this invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
[0019] The spectrophotometer of the disclosed embodiment is a 
spectrophotometer especially suitable for being mounted at one side of the printed 
1 5 sheets output path of a color printer to optically evaluate color imprinted output sheets 
as they move past the spectrophotometer, variably spaced therefrom, without having 
to contact the sheets or interfere with the normal movement of the sheets. In 
particular, it may be used to measure a number of color test patch samples printed by 
the printer on actual printed sheet output of the printer during regular or selected 
20 printer operation intervals (between normal printing runs or print jobs). These color 
test sheet printing intervals may be at regular timed intervals, and/or at each machine 
"cycle-up," or as otherwise directed by the system software. The spectrophotometer 
may be mounted at one side of the paper path of the machine, or, if it is desired to use 
duplex color test sheets, two spectrophotometers may be mounted on opposite sides of 
25 the paper path. 

[0020] Relatively frequent color calibration of a color printer is highly 
desirable, since the colors actually printed on the output media (as compared to the 
colors intended to be printed) can significantly change, or drift out of calibration over 
time, for various known reasons. For example, changes in the selected or loaded print 
30 media, such as differences paper or plastic sheet types, materials, weights, 

calendaring, coating, humidity, etc., or changes in the printer's ambient conditions, 
changes in the image developer materials, aging or wear of printer components, 
varying interactions of different colors being printed, etc. Printing test color patches 
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on test sheets of the same print media under the same printing conditions during the 
same relative time periods as the color print job being color-controlled is thus very 
desirable. 

[0021] It is thus also advantageous to provide dual -mode color test sheets, in 
5 which multiple color patches of different colors are printed on otherwise blank areas 
of each, or selected, banner, cover, or other inter-document or print job separator 
sheets. Different sets of colors may be printed on different banner or other test sheets. 
This dual use of such sheets saves both print paper and printer utilization time, and 
also provides frequent color calibration opportunities where the printing system is one 
10 in which banner sheets are being printed at frequent intervals anyway. 

[0022] An additional feature which can be provided is to tailor or set the 
particular colors or combinations of the test patches on a particular banner or other 
test sheet to those colors which are about to be printed on the specific document for 
that banner sheet, i.e., the document pages which are to be printed immediately 
1 5 subsequent to that banner sheet (the print job identified by that banner sheet). This 
can provide a "real time" color correction for the color printer which is tailored to 
correct printing of the colors of the very next document to be printed. 

[0023] The preferred implementations of the systems and features disclosed 
herein may vary depending on the situation. Also, various of the disclosed features or 
20 components may be alternatively used for such functions as gray scale balancing, 

turning on more than one illumination source at once, such as oppositely positioned 
LEDs, etc. 

[0024] It will be appreciated that these test patch images and colors may be 
automatically sent to the printer imager from a stored data file specifically designed 

25 for printing the dual mode banner sheet or other color test sheet page, and/or they may 
be embedded inside the customer job containing the banner page. That is, the latter 
may be directly electronically associated with the electronic document to be printed, 
and/or generated or transmitted by the document author or sender. Because the 
printed test sheet color patches colors and their printing sequence is known (and 

30 stored) information, the on-line spectrophotometer measurement data therefrom can 
be automatically coordinated and compared. 

[0025] After the spectrophotometer or other color sensor reads the colors of 
the test patches, the measured color signals may be automatically processed inside the 
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system controller or the printer controller to produce or modify the tone reproduction 
curve, as explained in the cited references. The color test patches on the next test 
sheet may then be printed with that new tone reproduction curve. This process may be 
repeated so as to generate further corrected tone reproduction curves. If the printer's 
5 color image printing components and materials are relatively stable, with only 

relatively slow long term drift, and there is not a print media or other abrupt change, 
the tone reproduction curve produced using this closed loop control system will be the 
correct curve for achieving consistent colors for at least one or even a substantial 
number of customer print jobs printed thereafter, and only relatively infrequent and 

1 0 few color test sheets, such as the normal banner sheets, need be printed. 

[0026] In addition to use in printers, it should be noted that color 
measurements, and/or the use of color measurements for various quality or 
consistency control functions, are also important for many other different technologies 
and applications, such as in the production of textiles, wallpaper, plastics, paint, inks, 

1 5 food products, etc. and in the measurement or detection of various properties of 

various materials, objects or substances. Thus, the invention may have applications in 
various such other fields where these materials, objects or substances are to be color 
tested, including both (1) applications in which color measurements are taken and 
applied in a closed loop control system and (2) applications in which the measurement 

20 result is not fed back into a control loop, but is used to generate a one-time output. 

[0027] Fig. 1 is a functional block diagram illustrating an exemplary 
embodiment of a coloring system 100 according to this invention. The coloring 
system 100 is connected to an input device 200 via a link 210. The input device 200 
inputs various information needed to implement the operations of the coloring system 

25 100, as described in more detail below, and may include a mouse, a keyboard, a touch- 
screen input device, a voice recognition-based input device, and/or any other known 
or later developed device usable for inputting information. The coloring system 100 
optionally is connected to an image data source 300 via a link 310. The connection to 
the image data source 300 is "optional" because it is required only for certain 

30 embodiments of the coloring system 100. 

[0028] For example, when the coloring system 100 is a marking device, such 
as a printer, the image data source 300 is required. However, when the coloring 
system 100 is a system for performing a coloring operation that does not require 



image data, the image data source 300 is not required. An example of a coloring 
operation that may not require image data is an operation of making a colored food 
product such as cereal or the like. 

[0029] The image data source 300 can be a digital camera, a scanner, or a 
locally or remotely located computer, or any other known or later developed device 
that is capable of generating electronic image data. Similarly, the image data source 
300 can be any suitable device that stores and/or transmits electronic image data, such 
as a client or a server of a network. The image data source 300 can be integrated with 
the coloring system 100, as in a digital copier having an integrated scanner. 
Alternatively, the image data source 300 can be connected to the coloring system 100 
over a connection device, such as a modem, a local area network, a wide area 
network, an intranet, the Internet, any other distributed processing network, or any 
other known or later developed connection device. 

[0030] It should also be appreciated that, while the electronic image data can 
be generated at the time of printing an image from an original physical document, the 
electronic image data could have been generated at any time in the past. Moreover, 
the electronic image data need not have been generated from the original physical 
document, but could have been created from scratch electronically. The image data 
source 300 is thus any known or later developed device which is capable of supplying 
electronic image data over the link 310 to the coloring system 100. The link 310 can 
thus be any known or later developed system or device for transmitting the electronic 
image data from the image data source 300 to the coloring system 100. 

[0031] Further, it should be appreciated that the links 210 and 3 1 0 can be a 
wired, wireless or optical link to a network (not shown). The network can be a local 
area network, a wide area network, an intranet, the Internet, or any other distributed 
processing and storage network. 

[0032] The coloring system 100 includes a coloring device 120, a sensor array 
130, a color revision device 140, a memory 150, a controller 160 and a spectral curve 
determination system 170, which are interconnected by a data/control bus 190. The 
spectral curve determination system 170 includes a reference database 172 and a 
spectral curve output device 174. 

[0033] The coloring device 120 may be, for example, a print engine/printing 
head or marking engine/marking head, when the coloring system 100 is a printer or 
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other marking device. The coloring device 120 may be, for example, a colorant 
dispenser that dispenses a colorant onto an object or into a mixture. In short, the 
coloring device 1 20 may be any known or later developed device that directly or 
indirectly controls the final appearance of an object, material or substance. 
5 [0034] The sensor array 130 includes multiple illuminants, such as LEDs, 

lasers or the like, arranged around a central photodetector (not shown), or arranged in 
correspondence to a plurality of photodetectors or photosites as described in above- 
mentioned co-pending U.S. Application No. 09/862,247, U.S. Application No. 
09/863,042, and/or U.S. Application No. 09/888,791. The illuminants will be referred 

10 to hereafter as LEDs for convenience. The number of LEDs may be any number 

greater than three, when a single photosensor is used, or may be as low as two when 
multiple photosites or photosensors are used. A larger number of LEDs gives greater 
accuracy, but it costs more to include more LEDs, and thus there are practical limits to 
the number of LEDs included in the sensor array 130, especially since an object of this 

1 5 invention is to provide a low-cost spectrophotometer. Therefore, the number of LEDs 
is preferably from about 8 to about 16. 

[0035] Each LED is selected to have a narrow band response curve in the 
spectral space. Therefore, for example, ten LEDs would correspond to ten 
measurements in the reflectance curve. The LEDs, or other multiple illuminant based 

20 color sensor equivalent, e.g., lasers, are switched on one at a time as, for example, the 
measured media is passed through a transport of a printer. The reflected light is then 
detected by the photodetector and the corresponding voltage integrated and 
normalized with a white tile. The normalization may be performed periodically. For 
the normalization, use of a white tile calibration look-up table, which is stored in 

25 memory 1 50, is a standard practice in the color measurement industry. When the 
white tile calibration look-up table is used, the detector output is normalized to 
between 0 to 1 in accordance with the following equation: 

where V" is the black measurement sensing system offset of the 1 th LED, Vf* is the 
30 white tile measurements, V, is the sensor detector output, and i?, w is the reflectance 
spectra of the white tile at the mean wavelength of the i th LED. Any other known or 
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later developed method for normalization may alternatively be used. F m may be 
compensated for temperature variation. 

[0036] The color revision device 140 calibrates the output of the coloring 
device 120 in accordance with information obtained from the spectral curve output 
5 device 174 of the spectral curve determination system 170. This calibration may be 
performed as often as necessary or desired to maintain a desirable output of the 
coloring device 120. 

[0037] The memory 150 may serve as a buffer for information coming into or 
going out of the coloring system 100, may store any necessary programs and/or data 

10 for implementing the functions of the coloring system 100, and/or may store data at 
various stages of processing. The above-mentioned white tile lookup table may be 
stored in the memory 1 50 if desired. The reference database 172, described in more 
detail below, may also be stored in the memory 150 if desired. Furthermore, it should 
be appreciated that the memory 150, while depicted as a single entity, may actually be 

1 5 distributed. Alterable portions of the memory 150 are, in various exemplary 

embodiments, implemented using static or dynamic RAM. However, the memory 150 
can also be implemented using a floppy disk and disk drive, a writeable optical disk and 
disk drive, a hard drive, flash memory or the like. The generally static portions of the 
memory 1 50 are, in various exemplary embodiments, implemented using ROM. 

20 However, the static portions can also be implemented using other non-volatile memory, 
such as PROM, EPROM, EEPROM, an optical ROM disk, such as a CD-ROM or 
DVD-ROM, and disk drive, flash memory or other alterable memory, as indicated 
above, or the like. 

[0038] The controller 160 controls the operation of other components of the 
25 coloring system 100, performs any necessary calculations and executes any necessary 
programs for implementing the processes of the coloring system 100 and its individual 
components, and controls the flow of data between other components of the coloring 
system 100 as needed. 

[0039] The spectral curve determination system 170 determines and outputs 
30 spectral curves. Specifically, the spectral curve output device 174 outputs spectral 
curves based on a plurality of spectra which are determined by the controller 160 
based on information from the reference database 172, described in more detail below, 
and the output of the sensor array 130 from different color targets. 
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[0040] To obtain an output similar to that of a reference spectrophotometer, 
such as a Gretag spectrophotometer, it is necessary to convert the readings from the 
sensory array 130 to reflectance spectra. A Gretag spectrophotometer outputs 36 
spectral reflectance values, evenly spaced at lOnm over the visible spectrum (e.g., 380 
5 nm to 730 nm). Therefore, in the following examples, the readings from the sensor 
array 130 are converted to 36 reflectance values. In other words, when there are 10 
LEDs in the sensor array 130, the LEDs are sequentially switched, readings (typically 
voltage readings) are collected from the photodetector for each respective LED, and 
the 10 readings (voltages) from the sensor array 130 are converted to 36 reflectance 

1 0 values per color. If a multiple photosite sensor is used, it will be appreciated that a 
desired number of outputs, for example 1 0 outputs, will be obtained from smaller 
number of LEDs, for example 3 or 4 LEDs. 

[0041] The reference database 172 is generated by measuring the reflectance 
spectra of some set of reference colors, with an accurate reference spectrophotomer, 

1 5 such as a Gretag spectrophotometer, and their corresponding LED sensor outputs, 

with the sensor array 130. In general, the more densely populated the database is, i.e., 
the more reference colors used, the better the resulting accuracy. In one exemplary 
reference database, about 2000 colors were used. Furthermore, even spacing of the 
reference colors in the color space gives greater accuracy. The data stored in the 

20 reference database 1 72 will be referred to hereafter as the training samples. 

[0042] It should be understood that each of the circuits shown in Fig. 1 can be 
implemented as portions of a suitably programmed general purpose computer. 
Alternatively, each of the circuits shown in Fig. 1 can be implemented as physically 
distinct hardware circuits within an ASIC, or using a FPGA, a PDL, a PLA or a PAL, 

25 or using discrete logic elements or discrete circuit elements. The particular form each 
of the circuits shown in Fig. 1 will take is a design choice and will be obvious and 
predictable to those skilled in the art. 

[0043] An exemplary algorithm that may be implemented by the controller 
160 for determining spectra based on the reference database 172 and the output of the 

30 sensor array 1 30 is described below. The following algorithm is a dynamic, Least 
Squares-based spectral reconstruction algorithm, hereafter referred to as a DLS 
algorithm. One characteristic of the DLS algorithm is that, in the determination of 
spectra, greater importance is given to the data from the training samples in the 
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neighborhood of the color sample under measurement. This is done using linear 
operators and basis vectors. 

[0044] In the following description, the number of LEDs included in the 
sensor array 1 30 is assumed to be 1 0. Those skilled in the art will appreciate how to 
apply the algorithm to sensor arrays with more or fewer LEDs. 

[0045] Furthermore, it should be appreciated in this context that, in general, 
algorithms applicable to generation of a tone reproduction curve are not applicable to 
generation of a spectral curve. One reason for this is that, while the first and last 
values in a tone reproduction curve are known (i.e., they are [0,0] and [255, 255]), the 
same cannot be said of spectral curves generated using LED sensors, because the 
LEDs at the opposite ends of the spectrum (i.e., the blue and red LEDs) are not 
monochromatic. 

[0046] First, the sensor transfer function, i.e., the information included in the 
reference database 172, is a mapping from reflectance spectra (as measured by a 
reference spectrophotometer) to sensor outputs (as measured by the sensor array 130) 
formed by a set of N spectra to voltage measurements, denoted as 

Q = [S, S 2 ... S N }e R"* N > Z = [Z, Z 2 ... ZJ e R'* N ...(2) 

where S x S 2 ...S N are the vector elements containing the N spectral curves, each 
curve containing 36 elements, i.e., reflectance values {n = 36), and Z, Z 2 ... Z N are 
the vector elements from the LED sensor outputs (in volts), each having ten 
components (/ = 10) when a 10-LED spectrophotomer is used. Here, each curve 
contains 36 elements because a Gretag spectrophotometer, which outputs 36 values, is 
used. If a different spectrophotometer is used which has a different number of 
outputs, n will be a correspondingly different number. Z, Z 2 ... Z N are each a vector 
including 10 normalized voltages corresponding to the 10 LED color sensor outputs 
for a given color. R indicates the set of real numbers. JVis a predetermined number 
based on certain color gamut specifications for the color sensor array 130. Generally, 
the larger the gamut, the larger will be N . As an example, N may be about 2000. 

[0047] The value of / discussed above depends on the number of sensor 
outputs. It will be appreciated that when a multiple photosite sensor is used, / will not 
be equal to the number of LEDs. 
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[0048] Given V m , which represents the set of 10 LED sensor voltages 
corresponding to an unknown (target) color with true spectrum S , a goal of the 
algorithm is to obtain an estimate of the unknown, true spectrum S . The relationship 
between color spectrum and LED sensor outputs for a given color may be assumed to 
be linear affine, quadratic affme or cubic affine etc. Therefore, the general expression 
for S may be defined as 

S = AV ...(3) 
where V — [V m 1] is an 1 1 x 1 vector for a linear affine model with 10 LEDs in the 
sensor array 130. V is a vector formed by augmenting the normalized V m voltage 
vector with scalar value 1 to include an affine term. Thus, for example, if 12 LEDs 
were used in the sensor array 130, V would be a 13 x 1 vector for a linear affine 
model. A is an unknown conversion matrix of size 36 x 1 1 for a 10 LED sensor and a 
linear affine model. Those skilled in the art will appreciate that, if quadratic or other 
terms are included, then the number of elements in vector V and matrix A 
correspondingly increase to greater than 1 lxl and 36x1 1 respectively. For example, if 
a quadratic term is included, then the number of elements in vector V will be 66 x 1 
and the number of elements in matrix A will be 36 x 66. 

[0049] To compute matrix A optimally for obtaining an estimated spectrum 
S for a color measured by sensor array 130, whose output is the augmented vector 
v = W m 1] , a weighted least square criteria is used. The optimal solution for 
matrix A , called A" , is obtained by minimizing the objective function defined as 

A* = arg min J = arg min J] w(i)\\S, - AZ, f ...(4) 

where 

wii)= lW7i ' ^H^-z.I -.(5) 

p is an integer number greater than or equal to 2, and s is a small positive constant. 
/ represents the index for training samples from the reference database (S, ,Z, ). 
Appropriate values for p and e may be empirically determined. Examples of 
appropriate values chosen for an 8- or 10-LED sensor are p = 4 and e = 1 x 10" 5 . 
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The solution to the above optimization problem can be easily obtained by setting the 
gradient of J with respect to A equal to zero. This results in 
A*=QP~ X ...(6) 

where 

N 

5 2 = ][>(/)S,Z, 7 ' ....(7) and 

P = f i w(i)Z l Zf ...(8) 

/=i 

[0050] Once A* is computed, the estimated spectrum S for the target color 
measured by the sensor array 130 can be obtained by 

10 

S = A'V ...(9) 

[0051] The foregoing algorithm determines one spectrum, obtained from one 
target color under measurement. If spectral curves are to be generated for a plurality of 

1 5 target colors, then matrix A* is computed for the plurality of target colors, using 
equations (5) - (9) for each color with corresponding V vector. 

[0052] Fig. 2 is a flowchart illustrating a first exemplary method of 
determining spectra according to this invention. Beginning in step SI 000, the process 
continues to step SI 100, where training samples are entered from a reference 

20 database, and appropriate values are entered for p , s and N . The process then 

continues to step SI 200 and receives a sensor reading, such as a sensor voltage, from 
each illuminant in a sensor array. Continuing to step SI 300, the process determines 
V by normalizing the sensor readings obtained from the sensor array, e.g., based on a 
white tile calibration look-up table. 

25 [0053] Next, in step S 1 400, the process determines matrix A* in accordance 

with the algorithm described above. The process then continues to step S 1 500, where 
a spectrum is determined as S = A 'V . Continuing to step SI 600, it is determined 
whether all color samples have been measured. If not all the color samples have been 
measured, the process continues to step SI 700. Otherwise, the process jumps to step 

30 S1800. 

[0054] In step SI 700, the next color sample is selected. Steps SI 200 - SI 600 
are then repeated. When all color samples have been measured, the process goes to 
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step SI 800 and outputs the full reflectance spectra, i.e., the spectral curve, of the color 
samples. Finally, the process ends in step SI 900. 

[0055] One disadvantage of the above-described algorithm is computational 
complexity due to the need for a matrix inversion, seen in the step of determining 
A* = QP~ l , every time a new spectrum is required for each color measured by the 
sensor array 130. Matrix inversion is time consuming and can even be an ill posed 
problem, since P and P must be positive definite matrices. A modified DLS 
algorithm which overcomes this problem is described below. 

[0056] Fig. 3 is a flowchart illustrating a second exemplary method of 
determining spectra according to this invention. Beginning in step S2000, the process 
continues to step S2050, where N training samples are entered from a reference 
database, and appropriate values are entered for p , s and N . Additionally, B(0) is 
set at a 2 I, where a » 1 and / is the identity matrix of size 1 1 x 1 1 for a 10-LED 
sensor, D(0) is a zero matrix of size 36 x 1 1 , and i, the index of recursive iteration, is 
set at 0. An appropriate value for a may be determined empirically. An example of 
an appropriate value for a is 1 x 10 5 . 

[0057] The process then continues to step S2100 and receives a sensor 
reading, such as a sensor voltage, from each illuminant in a sensor array. Continuing 
to step S21 50, the process determines V by normalizing the sensor readings obtained 
from the sensor array, e.g., based on a white tile calibration look-up table. The 
process then continues to step S2200. 

[0058] In step S2200, the process sets i = i + 1 , and then continues to step 
S2250. In step S2250, the process computes B(i + 1) and D(i + 1) as shown, where 
B(i) is an 1 1 x 1 1 matrix and D(i) is a 36 x 1 1 matrix, for a 10 LED sensor. The 
process then continues to step S2300 and determines whether i = N . If / = N , the 
process continues to step S2350. Otherwise, the process returns to step S2200. Thus, 
the process performs a recursive loop until i = N . 

[0059] In step S2350, an estimated spectrum S is determined as 
S = D(N)B(N)V . Continuing to step S2400, it is determined whether all color 
samples have been measured. If not all the color samples have been measured, the 
process continues to step S2450. Otherwise, the process jumps to step S2500. 

[0060] In step S2450, the next color sample is selected. Steps S2100 - S2400 
are then repeated. When all color samples have been measured, the process goes to 
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step S2500 and outputs the full reflectance spectra, i.e., the spectral curve, of the color 
samples. Finally, the process ends in step S2550. 

[0061] As seen above, the method of Fig. 3 does not require a matrix 
inversion. However, the method of Fig. 3 does require a recursive loop, which 
5 increases the processing time. Thus, for some applications, the method of Fig. 2 may 
be preferred over the method of Fig. 3, while for other applications, the method of Fig. 
3 may be preferred over the method of Fig. 2. Those skilled in the art will be able to 
determine which method is more suitable for a given application. 

[0062] Fig. 4 is a functional block diagram illustrating an exemplary 

1 0 embodiment of a color detection system 500 according to this invention. The color 
detection system 500 includes an input/output interface 1 10, a sensor array 130, a 
controller 150, a memory 160 and a reference database 172, which may be identical to 
the corresponding elements of Fig. 1, interconnected by a data/control bus 590. The 
color detection system 500 is connected to a user input device 200 via a link 210, 

15 similar to the input device 200 and link 210 described above in conjunction with Fig. 
1 . The color detection system 500 is also connected to a data sink 400 via a link 410 
which, like the links 2 1 0 and 3 1 0, can be a wired, wireless or optical link to a network 
(not shown). In general, the data sink 400 can be any device that is capable of 
outputting or storing the processed data generated by the color detection system, such 

20 as a printer, a copier or other image forming devices, a facsimile device, a display 
device, a memory, or the like. 

[0063] The color detection system 500 may be, or be included in, a portable or 
stationary unit designed specifically to measure color of a target object. In use, the 
color detection system 500 is positioned with the sensor array 130 facing the target 

25 object, the sensor array 130 is activated as described above, and then the above- 
described algorithm is executed by the controller 150, using data from the sensor array 
130 and the reference database 172, to obtain an estimated spectrum S of the target 
object. The estimated spectrum S is then output to the data sink 400. 

[0064] From the foregoing descriptions, it can be appreciated that, in 

30 embodiments, the invention may provide a calibration tool for scanners, printers, 

digital photocopiers, etc., and that, in embodiments, the invention may provide a color 
measurement tool designed to provide one-time color measurements of target objects. 
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[0065] The coloring system 100 of Fig. 1 and the color detection system 500 
of Fig. 4 are preferably implemented either on a single program general purpose 
computer or separate programmed general purpose computer, with an associated 
sensor array 130 (and coloring device 120, in the case of Fig. 1). However, the 
coloring system 100 and color detection system 500 can also be implemented on a 
special purpose computer, a programmed micro-processor or micro-controller and 
peripheral integrated circuit element, an ASIC or other integrated circuit, a digital 
signal processor, a hard-wired electronic or logic circuit such as a discrete element 
circuit, a programmable logic device such as a PLD, PLA, FPGA, PAL, or the like. In 
general, any device capable of implementing a finite state machine that is in turn 
capable of implementing the flowcharts shown in Fig. 2-3, or appropriate portions 
thereof, can be used to implement the spectral curve reconstruction device according 
to this invention. 

[0066] Furthermore, the disclosed methods may be readily implemented in 
software using object or object-oriented software development environments that 
provide portable source code that can be used on a variety of computer or workstation 
hardware platforms. Alternatively, appropriate portions of the disclosed coloring 
system 100 and the color detection system 500 may be implemented partially or fully 
in hardware using standard logic circuits or a VLSI design. Whether software or 
hardware is used to implement the systems in accordance with this invention is 
dependent on the speed and/or efficiency requirements of the system, the particular 
function, and the particular software or hardware systems or microprocessor or 
microcomputer systems being utilized. The processing systems and methods 
described above, however, can be readily implemented in hardware or software using 
any known or later developed systems or structures, devices and/or software by those 
skilled in the applicable art without undue experimentation from the functional 
description provided herein together with a general knowledge of the computer arts. 

[0067] Moreover, the disclosed methods may be readily implemented as 
software executed on a programmed general purpose computer, a special purpose 
computer, a micro-processor, or the like. In this case, the methods and systems of this 
invention can be implemented as a routine embedded on a personal computer or as a 
resource residing on a server or workstation, such as a routine embedded in a 
photocopier, a color photocopier, a printer driver, a scanner, or the like. The systems 
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and methods can also be implemented by physical incorporation into a software and/or 
hardware system, such as the hardware and software system of a photocopier or a 
dedicated image processing system. 

[0068] While the invention has been described in conjunction with the 
5 specific embodiments described above, many equivalent alternatives, modifications 
and variations may become apparent to those skilled in the art when given this 
disclosure. Accordingly, the exemplary embodiments of the invention as set forth 
above are considered to be illustrative and not limiting. Various changes to the 
described embodiments may be made without departing from the spirit and scope of 
10 the invention. 



